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Substantial anisotropies should occur in the hot expanding QCD plasma produced in 
relativistic heavy-ion collisions due to non-vanishing shear viscosity. We discuss the form of 
the real-time, hard thermal loop resummed propagator for static gluons in the presence of 
such anisotropies and the consequences for quarkonium binding. It has been predicted that 
the propagator develops an imaginary part due to Landau damping at high temperature. 
This should generate a much larger width of quarkonium states than the Appelquist-Politzer 
vacuum estimate corresponding to decay into three gluons. We argue that this might be 
observable in heavy-ion collisions as a suppression of the T(1S) —¥ e + e~ process. Lastly, we 
consider the heavy quark (singlet) free energy just above the deconfinement temperature. In 
the "semi-QGP", Fqq(R) at distances beyond 1/T is expected to be suppressed by 1/N as 
compared to an ideal plasma. 
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§1. The anisotropic QGP 



This section deals with kinetic non-equilibrium effects in an expanding QCD 
plasma. During the first few fm/c when the temperature is high it is the longitu- 
dinal Bjorken expansion^ 1 which matters most; for semi-peripheral collisions and/or 
close to the periphery of the fireball the transverse expansion may also be important 
but is neglected here for simplicity. As a consequence of the expansion, the particle 
momentum distribution in the local rest frame is anisotropic, if the scattering rate 
is finite: there is a net loss of particles with large \p z \ from the fluid cell and redistri- 
bution of the momenta requires time. We assume that the momentum distribution 
can be parameterized as follows P 



Hp) = /iso(^ P 2 + £(p-n) 2 ) ~ f iso (p) 
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(1-1) 



The second expression is valid to first order in the anisotropy parameter £. fi BO (p) 
is either a Bose or a Fermi distribution, respectively, and n = e z for longitudinal 
expansion. Such a correction 5f to the equilibrium distribution is also commonly 
employed in viscous hydrodynamics and has been argued to be consistent with az- 
imuthal flow coefficients observed at RHIC 3 ' (a transverse unit vector n = ex is 
used for calculating pj- distributions). Near equilibrium one can relate moments of 
5f, which are proportional to £, to the shearP 

Fig. [1] shows the time evolution of the ratio of the longitudinal to the trans- 
verse pressure in the central region of a high-energy heavy-ion collision (see, also, 
ref.ffl). This was obtained in refP 1 from a Boltzmann equation in relaxation time 



typeset using P7PTeX.cIs (Ver.0.9) 



2 



Adrian Dumitru 



1 

0.8 



0.6 
0.4 



0.2 







4o 


- 1 

= 10 








100 






















y / 


s' 

s' 










✓ 



















0.2 0.4 



0.8 1.6 3.2 
X (fm/c) 



6.4 



Fig. 1. Pressure anisotropy during the high-temperature stage of a heavy- ion collision^ for a low- 
viscosity plasma with rj/s = 0.1. The curves correspond to different initial anisotropics at 
T = 0.2 fm/c. 



approximation to all orders in the anisotropy parameter £. Near equilibrium, the 
assumed equilibration rate translates into a shear viscosity to entropy density ratio 
of rj/s = 0.1. However, even for such "strong coupling" conditions large pressure 
anisotropies are caused by the rapid longitudinal expansion during the early, high- 
temperature stage of a heavy-ion collision (r < 3 fm/c). Note that pi < Pt does 
not suppress transverse hydrodynamic flow^ but may reflect in the system-size and 
energy dependence of the transverse energy dE^jdy in the final stateP 

1.1. Real part of the quarkonium potential 

Anisotropies of the momentum distributions affect various processes in the plasma. 
In particular, rotational symmetry in the local rest frame is broken and the screen- 
ing length acquires an angular dependence. At high temperature and in the weak- 
coupling approximation it can be obtained explicitly from the "hard thermal loop" 
(HTL) resummed propagator for static electric gluons. To linear order in 

Re ^ 00 (p) = - T -U r (l - tml 2 *-f' n)2 2 /p2 ) . (1-2) 
jr + m D \ p + m D J 



Here, mo = gTyNj2> denotes the standard Debye mass of the equilibrium plasma. 
The one-gluon exchange potential follows essentially from the Fourier transform, 



Re V(r) = V iso (r) 1 + £ 



»2 "1 
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(1-3) 



Here, f = rmrj, cos 9 = r • n/r, and V; so (r) = — " s F exp(— f) is the well-known 
Debye-screened Coulomb potential. The potential (|T3p is valid for distances f < 1. 
At fixed T the ^-dependent correction in eq. (|l-3p reduces thermal screening effects 
as compared to an ideal plasma in local equilibrium (£ = 0). On the other hand, one 
can absorb this £ dependence to a large extent by a redefinition of the hard scale 
T(£) and/or of the Debye mass ttid(£).® The most useful approach with respect to 
applications to heavy- ion collisions would probably amount to adjusting the initial Tq 
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in such a way that the entropy (per unit of rapidity) in the final state remains fixed 
as £o is varied. The resulting implicit dependence To(£o) can be determined only 
from solutions of viscous hydrodynamics or transport theory (with gluon-number 
changing processes^). 
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Fig. 2. Left: Spectral function of pseudo-scalar bottomonium obtained from a potential modePJ 
with Voo(T) ~ 1/T. In this calculation, the width of the peaks has been put in by hand for 
numerical stability. Note the logarithmic scale. Right: Binding energy of P-state bottomonium 
(L — 1) as a function of temperature for an isotropic (£ = 0) and a moderately anisotropic 
(£ = 1) QGPp) The solid increasing line shows T itself; note that T c = 192 MeV was assumed 
here to normalize the abscissa. 



For semi-realistic estimates of the binding energies of charmonium and bot- 
tomonium states one needs to add the linear confining potential ~ o~r, where a ~ 
1 GeV/fm is the SU(3) string tension. Its temperature dependence to be used in the 
real-time formalism must presently be modelled. In the phenomenologically relevant 
range T/T c = 1 — 3 the "interaction measure" (e — 3p)/T 4 in QCD is large. The 
potential at infinite separation is thus sometimes modelled as^ 1 V 0O (T) ~ 2a /T with 
a ~ 0.08 GeV 2 a constant of dimension two. Such a model has in fact been proposed 
long ago in refP^* 



V(r) 



r mo 



1) 



or 



(1-4) 



which also provides a smooth interpolation to short distances. In this model the tem- 
perature dependence of the binding energy £"bind = H — — Iraq of small 

bound states such as the T actually turns out to originate mostly from ^(T).^ 
This is confirmed by the spectral function of pseudo-scalar bottomonium published 
in refP^ 1 which has been replotted on a logarithmic scale in fig. [2j while the ground 
state peak shows little temperature dependence, the continuum threshold decreases 
rapidly as T increases. Similar spectral functions have been shown in refP^ 

The anisotropy or viscosity affects the excited states more strongly than the 
compact 15 bottomonium state. In this model the binding energy at T/T c = 1.1, 
for example, increases by about 50%. Alternatively, if a "dissolution" temperature 
is defined from \Eb\ = T^ s then this increases from ~ 1.15T C when £ = to about 
1.35T C when £ = 1. Thus, excited states should melt less easily if the QGP exhibits 
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an anisotropy during the early stages of the expansion. Lastly, there is a splitting 
of the states with L • n = and L • n = ±1, respectively, since rotational symmetry 
is broken. Understanding the properties of excited states of the T is important for 
phenomenology because they contribute through feed down to the yield of the 15 
stated 

1.2. Imaginary part of the quarkonium potential 

At finite temperature, the quarkonium potential also acquires an imaginary 
partP^ at order g 2 Cp due to Landau damping of the exchanged nearly static gluon. 
It can again be obtained from the Fourier transform of the HTL-resummed real- 
time propagator (the "physical" component of the Schwinger-Keldysh representa- 
tion) for static Aq fields. Taking the imaginary part corresponds to cutting open 
one of the hard thermal loops of the HTL propagator and can be viewed, micro- 
scopically, as the scattering of the space-like exchanged gluon off a thermal gluonP^ 
g + (QQ) — > g + Q + Q- To order £ and in the leading log 1/f approximation ImV 
is given by^* 

Im V(r) = - 9 — r 2 log - - ) ■ (1-5) 

This corresponds to a decay width for a Coulomb ground state of 

167rTm 2 D f A g 2 C F M Q 

For £ = one can in fact evaluate the matrix element of Im V between Coulomb 
wave functions without resorting to the log 1/f >1 approximation, 

T r^ l-(2-K 2 ) 2 + 41ogi K _^_H}D h7) 
a s C F M 2 (1 - k 2 ) 3 ' a s C F Mq { ' 




Fig. 3. Left: Imaginary part of the real-time propagator due to Landau damping. Right: Thermal 
decay width as a function of temperature for an isotropic medium (£ = 0). 



The width obtained from (|1 T[) is shown in fig. [3j For temperatures accessible 
to the RHIC and LHC colliders, T r is on the order of 20 MeV - 100 MeV. This 
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can be compared to the T — > e + e~ decay width in vacuum which arises from the 
annihilation of the b, b quarks into di-electrons^]: r T ^ e + e - = 1.34 keVP^ Because 
the electromagnetic decay width is much smaller than the inverse lifetime of the 
QGP the actual T peak observed in the di-lepton invariant mass distribution would 
not be broadened as compared to vacuum. 

Nevertheless, T states which have been broken up into individual b and b quarks 
reduce the yield of di-leptons in the peak. One of the contributions to this process 
is 9 + (QQ) Q + Q dissociation^ by a thermal gluonPS'® However, it has been 
found that this process leads to rather small T dissociation rates. Indeed, refP^ 
argues that in the limit of loose binding, \Ef,\ <C T due to strong screening of the 
attractive interaction, it becomes more efficient to scatter a thermal gluon off one 
of the quasi-free 6-quarks thereby breaking up the bound state. The width (|1-7|) 
obtained from the imaginary part of the HTL propagator indicates that damping 
of the exchanged gluon in the heat bath also provides a large contribution to the 
thermal T — > b + b rate. 

At low temperature, once pions have formed, a non-zero (but exponentially 
small) width emerges due to 7r + T — > B + B corresponding to tunneling of the b, b 
quarks in the background field of the pionPU 1 

It is certainly interesting to compare to a strongly coupled theory. Using the 
gauge-gravity duality, the static potential (or Wilson loop)p3 and thermal effects at 
short distances 2 ^ 1 have been computed in M = 4 supersymmetric Yang- Mills at large 



(but finite^) t' Hooft coupling A = g 2 N and N 



oo. 



At T = 0, 



V QQ (r) 



4tt 2 V\ 
T(l/4) 4 r 



The ~ 1/r behavior follows from conformal invariance of the theory. Also, the 
potential is non-analytic in A. Clearly, the coupling should not be very large or else 
the properties of the resulting bound states are qualitatively different from the T 
etc. states of QCD (numerically, 4tt 2 / T '(1 / '4) 4 » 0.23). 

At finite temperature, the potential develops an imaginary part when the string 
dangling in the fifth dimension (with end points on our brane) approaches the black 
hole horizon which sets the temperature scale and "melts" P^ 1 Fluctuations near the 
"bottom" of the string should generate such an imaginary part of the Nambu-Goto 
action at even lower temperatures already,^ 



llm V, 



QQ\ 



ttVX b 
48^/2 «o 



45 



a T 



(1-9) 



where denotes the unperturbed Coulomb ground state wave function, ao = 
r(l/4) 4 /27r 2 \/Am Q is the "Bohr radius" for the Maldacena potential (fT8|) and 
b = IF (3/4) jypH r(l/4) ~ 0.38 is a numerical constant. Here the width decreases 
with the quark mass and with the t' Hooft coupling approximately as -Tqq ~ 1/A m,Q ; 



*' It is therefore proportional to the square of the quarkonium wave function at the origin, 
similar to the hadronic decay into three gluons discussed by Appelquist and Politzer^ but unlike 
the thermal QCD width discussed above. 
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it increases rapidly with the temperature, ~ T 4 . For uiq = 4.7 GeV, T = 0.3 GeV, 
y/X = 3 one finds T r ~ 50 MeV. The (QQ) -> (Qq) (Qq) breakup due to string 
splitting at has been considered in ref.® 

The suppression of the yield of di-leptons from Y(1S) decays in the final state 
should be significant.®'® Neglecting "regeneration" of bound states from b and 
b quarks in the medium, the number of T mesons in the plasma which have not 
decayed into unbound b and b quarks up to time r after the collision is 

N(t) ~ iV exp ^- J dtPr(t)j . (1-10) 

This solution assumes that r~f(T(t)) is a slowly varying function of time. The initial 
number of Y states may be estimated from the multiplicity in p+p collisions times 
the number of binary collisions at a given impact parameter: No ~ N co \\Np p . Thus, 
the "nuclear modification factor" Raa for the process Y — > £ + £~ is approximately 
given by Raa{Y — > — exp(— Ptt), where F denotes a suitable average of 

r(T) over the lifetime of the quark-gluon plasma. Due to the strong temperature 
dependence of the width, this average is dominated by the early stage. Experimental 
measurements of Raa(T — > at RHIC appear to indicate a suppression® but 

so far it has not been possible to disentangle various Y states. 



§2. The static gluon propagator in the "semi-QGP" 



This section is about the HTL propagator in Euclidean time for temperatures 
just above deconfinement. Hidaka and Pisarski suggested that for T not far above 
T c hard modes in the QCD plasma may still be weakly coupled but propagating in 
a non-perturbative background Aq field.® The expectation value of the Polyakov 
loop is parameterized as 

L = exp(-i^ , Q = gA , £ = -^tr L . (2-1) 

At high T the eigenvalues of the matrix L approach and so its normalized trace 




Fig. 4. Left: Illustration of "repulsion" of the eigenvalues of the Polyakov loop for N = 3 colors. 
Right: Singlet free energy versus distance R, data fromP^ 
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£ — > 1. On the other hand, as T — > T c , the eigenvalues "repell" and eventually, at 
T c — 0, distribute uniformly over a circle such that I = 0. Note that £(T c + 0) ~ 0.4 in 
SU(3) Yang-Mills, 30 ' far from unity, and hence that Aq ~ T/g is non-perturbatively 
large. It is in this sense that the "semi-QGP" just above T c is a weakly coupled but 
non-perturbative phase. 

Gluons which couple to the background field acquire an additional "mass"; the 
static propagator in the semi-QGP (summed over colors) becomes, schematicalhF^l 

i(, !><*) L(0» ~ tC, J ^ £ p2 + mi ; {Ql _ Qk? ■ (2-2) 

The -/V diagonal gluons are screened only over large distances of order 1/mo, where 
m 2 D = 0{g 2 T 2 N) is the usual Debye mass. The iV 2 — N off-diagonal gluons corre- 
sponding to I 7^ k, however, acquire a mass of order T when Aq — Aq ~ T/g. The 
heavy-quark free energy Fqq(R) /T defined from the correlator of Polyakov loops 
(which may not necessarily be identified with the potential in real time) is expected 
to show the following qualitative behavior (assuming, for simplicity, large N): 

RT < 1 

1 < RT < - 
9 

- < RT 
9 

In the intermediate region f|2 -4[) the iV 2 — ./V heavy off-diagonal gluons have decoupled 
and Fqq is suppressed by 1/N as compared to its behavior at short distances. In 
fig. H] we show the singlet free energy for SU(3) Yang- Mills just above T c obtained 
numerically on 32 3 x4 lattices.^ Indeed, it does appear to show a two-slope behavior 
with a transition at RT ~ 2. On the other hand, near T c the large distance behavior 
could be affected by finite volume artifacts^" 1 simulations on larger lattices may 
hopefully become available in the future. If the behavior seen in the present data is 
confirmed then this may indicate that just above T c the properties of large, excited 
quarkonium states may be modified as compared to predictions from a standard 
potential model such as (|l-4p . 
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